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Abstract. Grids of models for rotating stars are con- 
structed in the range of 9 to 120 M at solar metallic- 
ity. The following effects of rotation are included: shellu- 
lar rotation, new structure equations for non-conservative 
case, surface distorsions, increase of mass loss with rota- 
tion, meridional circulation and interaction with horizon- 
tal turbulence, shear instability and coupling with ther- 
mal effects, advection and diffusion of angular momentum 
treated in the non-stationary regime, transport and dif- 
fusion of the chemical elements. 

Globally we find that for massive stars the effects of 
rotation have an importance comparable to those of mass 
loss. Due to meridional circulation the internal rotation 
law fi(r) rapidly converges, in 1-2 % of the MS lifetime, to- 
wards a near equilibrium profile which then slowly evolves 
during the MS phase. The circulation shows two main 
cells. In the deep interior, circulation rises along the po- 
lar axis and goes down at the equator, while due to the 
Gratton-Opik term it is the inverse in outer layers. This 
external inverse circulation grows in depth as evolution 
proceeds. We emphasize that a stationary approximation 
and a diffusive treatment of meridional circulation would 
be unappropriate. After the MS phase, the effects of core 
contraction and envelope expansion dominate the evolu- 
tion of the angular momentum. 

The surface velocities decrease very much during the 
MS evolution of the most massive stars, due to their high 
mass loss, which also removes a lot of angular momen- 
tum. This produces some convergence of the velocities, 
but not necessarily towards the break-up velocities. How- 
ever, stars with masses below ~12 M with initially high 
rotation may easily reach the break-up velocities near the 
end of the MS phase, which may explain the occurence 
of Be-stars. Some other interesting properties of the rota- 
tional velocities are pointed out. 

For an average rotation, the tracks in the HR dia- 
gram are modified like a moderate overshoot would do. 
In general, an average rotation may increase the MS life- 
time up to about 30 %; for the helium-burning phase the 



effects are smaller and amount to at most 10 %. From 
plots of the isochrones, we find that rotation may increase 
the age estimate by about 25 % in general. However, for 
stars with M J> 40 M Q and fast rotation, a bluewards 
"homogeneous-like" track, with important He- and N- 
enrichments, may occur drastically affecting the age esti- 
mates for the youngest clusters. Rotation also introduces a 
large scatter in the mass-luminosity relation: at the same 
logged an d logTcfj, differences of masses by 30 % may eas- 
ily occur, thus explaining what still remains of the alleged 
mass discrepancy. 

Rotation also brings significant surface He- and N- 
enhancements, they are higher for higher masses and ro- 
tation. While it is not difficult to explain very fast ro- 
tators with He- and N-excesses, the present models also 
well account for the many OB stars exhibiting surface en- 
richments and moderate or low rotation, (cf. Herrero et 
al. 1992, 200C). These stars likely result from initially fast 
rotators, which experienced mixing and lost a lot of angu- 
lar momentum due to enhanced mass loss. The comparison 
of the N-excesses for B- and A-type supergiants supports 



the conclusion by Venn (1995a, 1999), that these enrich 



ments mostly result from mixing during the MS phase, 
which is also in agreement with the results of Lyubimkov 
( |1996|) . 
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1. Introduction 

Due to the many well known successes of the standard the- 
ory of stellar evolution, rotation has generally been con- 
sidered as a secondary effect. This was justified in many 
cases. However, since some years a number of very signi- 
ficant discrepancies between model predictions and obser- 
vations have been found for massive stars, and also for red 
giants of lower masses. T hese dif ficulties have been listed 
and examined (cf. Maeder 1995a). Recently, large excesses 
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of [N/H] have been found in A-type supergiants, particu- 
larly in the SMC (Venn et al. |1998| ; Venn |1999| ), where the 
excesses may reach up to an order of magnitude with re- 
spect to the average [N/H] local value in the SMC. These 
excesses, not predicted by current models, are the signa- 
tures of mixing effects over the entire star. This extensive 
mixing changes the size of the reservoir of nuclear fuel 
available for evolution, and thus the lifetimes, the tracks 
and all the model outputs (cf. also Langer 1992; Langer 
1997; Heger et al. 2000). Also the chemical abundances, 
the yields and the final stages may be modified. Rotation 
appears as a natural driver for this mixing or, at least, as 
one of the first mechanisms whose consequences for mix- 
ing has to be explored. This is especially true for massive 
stars which are known to be fast rotators (see e.g. Penny 
1996; Howarth et al. 1997). 

The physical effects of stellar rotation are numerous. 
The basic equations of stellar structure need to be mo- 
dified. The meridional circulation and its interaction with 
the horizontal turbulence, the diffusion effects produced 
by shear turbulence, the inhibition of the transport me- 
chanisms by the //-gradients, the transport of the angular 
momentum and of the chemical elements, the loss of mass 
and angular momentum at the surface, the enhancement 
of mass loss by rotation, etc. . . are some of the effects to 
be included in realistic models. In addition to the model 
physics, rotation also brings a number of new numerical 
problems in the stellar code, such as the 4 th order equation 
for the transport of the angular momentum, which has to 
be coupled to the equations of stellar structure. 

In this work, we apply the investments in the treatment 
of the physical and numerical effects of rotation made in 
the previous works and we explore their consequences for 
the outputs of stellar models. In Sect. 2, we describe the 
various effects considered. The non-rotating stellar models 
are briefly discussed in Sect. 3. The evolution of the inter- 
nal rotation during the evolution is examined in Sect. 4. 
Sect. 5 discusses the effects of rotation on the evolutionary 
tracks in the HR diagram, on the lifetimes and isochrones. 
The evolution of the rotational velocities at the stellar sur- 
face is discussed in Sect. 6 and the effects on the surface 
abundances are analysed in Sect. 7. 

2. Physical ingredients of the models 

Let us briefly summarize here the basic physical ingre- 
dients of the numerical models of rotating stars we are 
constructing here. 

2.1. Shellular rotation 

The differential rotation which results from the evolution 
and transport of the angular momentum as described by 
Eq. (3) below, makes the stellar interior highly turbulent. 
The turbulence is very anisotropic, with a much stronger 
geostrophic-like transport in the horizontal direction than 



in the vertical one (Zahn 1992), where stabilisation is 
favoured by the stable temperature gradient. This strong 
horizontal transport is characterized by a diffusion coef- 
ficient -Dh, which is quite large as will be shown below. 
The horizontal turbulent coupling favours an essentially 
constant angular velocity O on the isobars. This rotation 
law, constant on shells, applies to fast as well as to slow 
rotators. As an approximation, it is often represented by 
a law of the form £1 = Vl(r) (Zahn 1992; see also Endal 



and Sofia |1976[) 



2.2. Hydrostatic effects 

In a rotating star, the equations of stellar structure need to 



be modified (Kippenhahn and Thomas 197C). The usual 
spherical coordinates must be replaced by new coordinates 
characterizing the cquipotcntials. The classical method 
applies when the effective gravity can be derived from a 
potential ^ = $ — ^fl 2 r 2 sin 2 9, i.e. when the problem is 
conservative. There, $ is the gravitational potential which 
in the Roche approximation is $ = — r . If the rotation 
law is shellular, the problem is non-conservative. Most 
existing models of rotating stars apply, rather inconsis- 
tently, the classical scheme by Kippenhahn an d Th omas. 
However, as shown by Meynet and Maeder ( 1997 ), the 
equations of stellar structure can still be written consis- 
tently, in term of a coordinate referring to the mass in- 
side the isobaric surfaces. Thus, the problem of the stellar 
structure of a differentially rotating star with an angular 
velocity Q = f2(r) can be kept one-dimensional. 



2.3. Surface conditions 

The distribution of temperature at the surface of a rota- 



ting star is described by the von Zeipel theorem (1924). 
Usually, this theorem applies to the conservative case and 
states that the local radiative flux F is proportional to the 
local effective gravity g e s , which is the sum of the gravity 
and centrifugal force, 



HP) 



9eff 



(1) 



with M*(P) = M(l - a^7=); H p ) is the luminosity on an 
isobar and p the mean internal density. The local T e g on 
the surface of a rotating star varies like T e ff($) ~ g c ff($) 1 ■ 
We define the average stellar T e g by T c 4 ff = L/(aS(Cl)), 
where a is Stefan's constant and 5(f2) the total actual 
stellar surface. Of course, for different orientation angles 
i, the emergent luminosity, colours and spectrum will be 



different (Maeder and Peytremann 1970). In the case of 
non-conservative rotation law, the corrections to the von 
Zeipel theorem depend on the opacity law and on the de- 
gree of differential rotation, but they are small, i.e. < 1% 



in curre nt ca ses of shellular rotation (Kippenhahn 1977 



Maeder 1999 ). Wether or not the star is close to the Ed- 
dington limit, the von Zeipel theorem keeps the same form 
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as the one given by Maeder (1999, see also Maeder & 
Meynet 2000a). In this last work, it is in particular shown 
that the expression for the Eddington factor in a rotating 
star needs to be consistently written and that it contains 
a term depending on rotation. 



2. 5. Transport of the angular momentum 

For shellular rotation, the equation of transport of angu- 
lar momentum in the vertical direction is in lagrangian 
coordinates (cf. Zahn 1992 ; Maeder and Zahn 1998] ) 



2.4- Changes of the mass loss rates M with rotation 

Observationally, a growth of the mass flux of OB stars 
with rot ation , i.e. by 2-3 powers of 10, was fou nd by 
Vardya ( 1985 ), while Nieuwenhuijzen and de Jager ( 1988 ) 
concluded that the M-rates seem to increase only slightly 
with rotation for O- and B-type stars. On the theoreti- 
cal side, Friend and Abbott (1986) find an increase of the 
M-rates which can be fitted by the relation (Langer 1998 ) 



P 



dt 

l_d_ 

5r 2 dr 



1 d 



__ (pr 4 fif/(r)) + __( pI)r 4_ 



Or 



(3) 



H(r) is the mean angular velocity at level r. The vertical 
component u(r, 6) of the velocity of the meridional circu- 
lation at a distance r to the center and at a colatitude 9 
can be written 



M{v) = M(v = 0) 



1 - 



(2) u(r,9) = U(r)P 2 (cos 9), 



(4) 



with £ ~ 0.5 and i> cr it the equatorial velocity at break- 
up. This expression is used in most evolutionary mod- 
els and this is also what is done in the present work. 
The critical rotation velocity of a star is often written 
as v^ it — ^P^(l — r), where i? c b is the equatorial ra- 
dius at break-up velocity and T = L/Lgdd is the ratio 
of the stellar luminosity to the Eddington luminosity (cf. 
Langer 1997, 1998). Glatzel (1998) has shown that when 
the effect of gravity darkening is taken into account, the 
above expression for w cr ; t does not apply. Glatzel (1998) 
gives v^ lit — which we adopt here. The problem is 

now being further examined by Maeder & Meynet (2000a) , 
who critically discuss the Eddington factors, their depen- 
dence on rotation, the expression of the critical velocity, 
the dependence of the mass loss rates on rotation. These 
various new results will be applied in subsequent works, 
particularly for the study of the effects of rotation on the 
formation of W-R stars. Further improvements, based ei- 
ther on the observations or on the theory, to account for 
the anisotropic winds which selectively remove the angular 
momentum need also to be performed. 

The above expression gives the change of the mass loss 
rates due to rotation. As reference mass loss rates in the 
case of no rotation, we use the recent data by Lamers and 
Cassinelli ( |1996[ ) ; for the domain not covered by these au- 
thors we use the results by de Jager et al. ( ]1988 ). During 
the Wolf-Rayet phase we use the mass loss rates proposed 
by Nugis et al (1998) for the WNL stars (mean, dumping- 
corrected rates from radio data M (WNL) = 3 10~ 5 M Q 
y _1 ). For the WNE and WC stars we use the prescription 
devised by Langer (1989), modified according to Schmutz 
(1997) for taking into account the clumping effects in 
Wolf-Rayet stellar winds (M = 2.4 • 1O- 8 (M/M ) 2 - 5 M 
y _1 ). These mass loss rates are smaller by a factor 2-3 
than the mass loss rates used in our previous stellar grids 
(Schaller et al. 1992; Meynet et al. 1994). 



where P2(cos9) is the second Legendre polynomial. Only 
the radial term U(r) appears in Eq. (3). The quantity D 
is the total diffusion coefficient representing the various 
instabilities considered and which transport the angular 
momentum, namely convection, semiconvection and shear 
turbulence. As a matter of fact, a very large diffusion coef- 
ficient as in convective regions implies a rotation law which 
is not far from solid body rotation. In this work, we take 
D = -D s hcar in radiative zones, since as extra-convective 
mixing we consider shear mixing and meridional circula- 
tion. In case the outward transport of the angular momen- 
tum by the shear is compensated by an inward transport 
due to the meridional circulation, we obtain the local con- 
servation of the angular momentum. We call this solution 
the s tation ary solution. In this case, U(r) is given by (cf. 
Zahn |1992| ) 



U(r) 



5D on 
~n~dr~ 



(5) 



The full solution of Eq. (3) taking into account U (r) and 
D gives the non-stationary solution of the problem. In 
this case, f2(r) evolves as a result of the various transport 
processes, according to their appropriate timescales, and 
in turn differential rotation influences the various above 
processes. This produces a feedback and, thus, a self- 
consistent solution for the evolution of f2(r) has to be 
found. 

The transport of angular momentum by circulation has 
often been treated as a diffusion process (Endal and Sofia 
1976| Pinsonneault et al. |1989| Heger et al. [2000] ). From 



Eq. (3), we see that the term with U (advection) is functio- 
nally not the same as the term with D (diffusion) . Physi- 
cally advection and diffusion are quite different: diffusion 
brings a quantity from where there is a lot to other places 
where there is little. This is not necessarily the case for 
advection. A circulation with a positive value of U(r), i.e. 
rising along the polar axis and descending at the equator, 
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is as a matter of fact making an inward transport of an- 
gular momentum. Thus, we see that when this process is 
treated as a diffusion, like a function of ffi, even the sign 
of the effect may be wrong. 

The expression of U(r) given below (Eq. 12) involves 
derivatives up to the third order, thus Eq. (3) is of the 
fourth order, which makes the system very difficult to solve 
numerically. In practice, we have applied a Henyey scheme 
to make the calculations. Eq. (3) also implies four bound- 
ary c onditi ons. At the stellar surfac e, we take (cf. Talon 
et al. |l997t Denissenkov et al. |1999[) 



where D\ x is the coefficient of horizontal turbulence, for 
which the estimate is 



D h ~ \rU{r)\ 



(10) 



according to Zahn (1992). Eq. (8) expresses that the ver- 
tical advection of chemical elements is severely inhibited 
by the strong horizontal turbulence characterized by -Dh- 
Thus, the change of the mass fraction Xi of the chemical 
species i is simply 



dr 



and U(r) = 



and at the edge of the core we have 



an 

dr 



and fi(r) 



(6) 



(7) 



We assume that the mass lost by stellar winds is just em- 
barking its own angular momentum. This means that we 
ignore any possible magnetic coupling, as it occurs in low 
mass stars. This is not unreasonable in view of the nega- 
tive results about the detection of magnetic fields in mas- 
sive stars (Mathys |1999 ). It is interesting to mention here, 
that in case of no viscous, nor magnetic coupling at the 
stellar surface, i.e. with the boundary conditions (6), the 
integration of Eq. (3) gives for an external shell of mass 
AM (Maeder 1999, paper IV) 



Am4(^ 2 ) = -— pr 4 nU(r). 
at 5 



(8) 



This equation is valid provided the stellar winds are spher- 
ically symmetric (see paper IV), an assumption we do in 
this work. When the surface velocity approches the criti- 
cal velocity, it is likely that there are anisotropics of the 
mass loss rates (polar ejection or formation of an equato- 
rial ring) and thus the surface condition should be modi- 
fied according to the prescriptions of Maeder ( |1999| ) . For 
now, these effects are not included in these models. Their 
neglect should not affect too much the results presented 
here since the critical velocity is reached only in some rare 
circumstances. 

2. 6. Mixing and transport of the chemical elements 

A diffusion-advection equation like Eq. (3) should nor- 
mally be used to express the transport of chemical ele- 
ments. However, if the horizontal component of the tur- 
bulent diffusion -Dh is large, the vertical advection of the 
elements c an b e treated as a simple diffusion (Chaboyer 
and Zahn 1992 ) with a diffusion coefficient D c g. As em- 
phasized by Chaboyer and Zahn, this does not apply to 
the transport of the angular momentum. D c g is given by 

1 rU{r) | 2 



D eS 



30D h 



(9) 



dXi 
dt 



dM r 



(Wp) 2 D n 



dXj 
dM r 



dXj 
dt 



(11) 



The second term on the right accounts for composition 
changes due to nuclear reactions. The coefficient D mlx 
is the sum £> mix = Ashear + -Deff and D cS is given by 
Eq. (9). The characteristic time for the mixing of chcmi- 
cal elements is therefore i m i x ~ tJ— an d is not given by 
t C i lc ~ as has been generally considered (Schwarzschild 



1958). This makes the mixing of the chemical elements 
much slower, since D c g is very much reduced. In this con- 
text, we recall that several authors have reduced by large 
factors, up to 30 or 100, the coefficient for the transport 
of the chemical elements, with respect to the transport 
of the angular momentum, in order to better fit the ob- 
served surface compositions (cf. Heger et al. 2000). This 



reduction of the diffusion of the chemical elements is no 
longer necessary with the more appropriate expression of 
D D g given here. 

When the effects of the shear and of the meridional 
circulation compensate each other for the transport of the 
angular momentum (stationary solution, see Sect. 2.5), the 
value of U entering the expression for D c g is given by Eq. 
(5). 

2. 7. Meridional circulation 

Meridional circulation is an essential mixing mechanism 
in rotating stars and there is a considerable litterature 
on the subject (see ref. in Tassoul 199C). The velocity of 
the meridional circulation in the case of shellular rotation 



was derived by Zahn (1992). The effects of the vertical 
/i-gradient V M and of the horizontal turbulence on merid- 
ional circulation are very important and they were taken 



into account by Maeder and Zahn (1998). Contrarily to the 
conclusions of previous works (e.g. Mestel 1965; Kippen- 
hahn and Weigert 1990; Vauclair 1999| ), the /x-gradients 
were shown not to introduce a velocity threshold for the 
occurence of the meridional circulation, but to progres- 
sively reduce the circulation when V M increases. The ex- 



pression by Maeder and Zahn ( 1998 ) is 



G. Meynet, A. Maeder: Stellar Evolution with Rotation 



5 



U(r) 



pgCpT [V ad 
L 



(Ec. 



V + (tp/6)V»] 



(12) 



P is the pressure, Cp the specific heat, i?n and are 
terms depending on the SI— and /^-distributions respec- 
tively, up to the third order derivatives and on vari- 
ous th ermod ynamic quantities (see details in Maeder and 
Zahn, |1998|) . The term V M in Eq. (12) results from the 
vertical chemical gradient and from the coupling between 
the horizontal and vertical /x-gradients due to the horizon- 
tal turbulence. This term V M may be one or two orders 
of magnitude larger than V a d — V in some layers, so that 
U(r) may be reduced by the same ratio. This is one of the 
important differences introduced by the work by Maeder 
and Zahn ( 1998 ). Another difference is that the classical 
solution usually predicts an infinite velocity at the inter- 
face between a radiative and a semiconvective zone with 
an inverse circulation in the semiconvective zone. Expres- 
sion (12) gives a continuity of the solution with no change 
of sign from semiconvective to radiative regions. Finally, 
we recall that in a stationary situation, U(r) is given by 
Eq. (5), as seen above. 

2.8. Shear turbulence and mixing 

In a radiative zone, shear due to differential rotation is 
likely to be a most efficient mixing process. Indeed shear 
instability grows on a dynamical timescale that is of the 
order of the rotation period (Zahn 1992] ). The usual crite- 
rion for shear instability is the Richardson criterion, which 
compares the balance between the restoring force of the 
density gradient and the excess energy present in the dif- 
ferentially rotating layers, 



Ri 



(0.8836 Vl^) 2 



1 

<4' 



(13) 



where we have taken the average over an isobar, r is the 
radius and iV a( j the Brunt- Vaisala frequency given by 



gs_ 

Hp 



if 

-gVfi + V ad — V ra d 



(14) 



When thermal dissipation is significant, the restoring force 
of buoyancy is reduced and the instability occurs more 
easily, its timescale is however longer, being the thermal 
timescale. This case is referred to as "secular shear insta- 
bility" . The criterion for low Peclet numbers Pe (i.e. of 
large thermal dissipation, see below) has been considered 



by Zahn (1974), while the cases of ge neral P eclet numbers 
Pe have been considered by Maeder ( |l995b| ), Maeder and 
Meynet ( 1996 ), who give 



(0.8836 Vl^f^) 2 H P 



(V ad - V) 



r + i 



< 



(15) 



The quantity T = Pe/Q, where the Peclet number Pe is 
the ratio of the thermal cooling time to the dynamical 
time, i.e. Pe = ^ where v and £ are the characteristic 
velocity and length scales, and K = (4acT 3 ) / (3CpKp 2 ) is 
the thermal diffusivi ty. A discussion of shear-driven tur- 
bulence by Canuto ( 1998 ) suggests that the limiting Ri 
number may be larger than 4. 

To account for shear transport and diffusion in Eqs. (3) 
and (11), we need a diffusion coefficient. Amazingly, a 
great variety of coefficients Z? s hcar — have been de- 
rived and applied, for example: 



- 1. Endal and Sofia (1978) apply the Reynolds and 
the Richardson criterion by Zahn ( [1974 ). They estimate 
-Dshear from the product of the velocity scale height of the 
shear flow and of the turbulent velocities of cells at the 
edge of Reynolds critical number. 

- 2. Pinsonneault et al. (198£) notice that the amount 
of differential rotation permitted by the secular shears is 
proportional to a critical number, which they treat as an 
adjustable parameter. To account for the effects of the 
/x-gradient in mixing, of the loss of angular momentum, 
etc. . . they introduce several adjustable parameters in the 
equations for the transport of the chemical elements and 
of the angular momentum. 

- 3. Chaboyer et al. (1995a) use a coefficient derived 



from the velocity and path le ngth f rom Zahn (1974). Fol- 
lowing Pinsonneault et al. ( 1989| ), they also introduce 
two adjustable parameters for adjusting the transports of 
chemical elements and angular momentum respectively. 
We notice that thanks to the reduction of the diffusion 
produced by the horizontal turbulence (cf. Eq. 9 above), 
it is no longer necessary to arbitrarily reduce the vertical 
transport of the chemical elements. 

- 4. Zahn (1992) defines the diffusion coefficient cor- 



responding to the eddies which have the largest Pe num- 
ber so that the Richardson criterion is just marginally 
satisfied. However, the effects of the vertical /i-gradient 
are not accounted for and the expression only applies to 
low Peclet numbers. 

- 5. The same has been done by Maeder and Meynet 
( 1996p , who considers also the effect of the vertical \i- 
gradient, the case of general Peclet numbers and, in ad- 
dition they account for the coupling due to the fact that 
the shear also modifies the local thermal gradient. This 



coefficient has been used by Meynet and Maeder (1997) 
and by Denissenkov et al. ( |1999[ ). 

- 6. The comparisons of model results and observations 
of surface abundances have led many authors to conclude 
that the /x-gradients appear to inhibit the shear mixing 
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too much with respect to what is required by the observa- 



tions (Chaboyer et al. 1995a b; Meynet and Maeder 1997; 
Heger et al. [2000 ). Namely, the observations of O-type 
stars (Herrero et al. 1992, 1999) show much more He- 
and N-enrichments than predicted by the models which 
apply Richardon's criterion with the V M term. Thus, in- 
stead of using a gradient V M in the criterion for shear mix- 



ing, Chaboyer et al. (1995a) and Heger et al. (2000) write 
/^Vp with a factor = 0.05 or even smaller. This proce- 
dure is not satisfactory since it only accounts for a small 
fraction of the existing /^-gradients in stars. The problem 
is that the models depend at least as much (if not more) 
on than on rotation, i.e. a change of / M in the allowed 
range (between and 1) produces as important effects as 
a change of the initial rotational velocity. This situation 
has led to two other more physical approaches discussed 
below. Also Heger et al. (2000) introduce another factor 
f c to adjust the ratio of the transport of the angular mo- 
mentum and of the chemical elements like Pinsonneault 
et al. (1989). 



- 7. Talon and Zahn ( 1997 ) account for the hori- 
zontal turbulence, which has a coefficient and which 
weakens the restoring force of the gradient of /j, in the 
usual Richardson criterion. This allows some mixing of 
the chemical elements to occur as required by the obser- 



vations (Talon et al. 1997). 

- 8. Indeed, around the convective core in the region 
where the /i-gradient inhibits mixing, there is anyway 
some turbulence due to both the horizontal turbulence 
and to the semiconvective instability, which is generally 
present in massive stars. This situation has led to the 



hypothesis (Maeder 1997) that the excess energy in the 
shear, or a fraction a of it of the order of unity, is degraded 
by turbulence on the local thermal timescale. This pro- 
gressively changes the entropy gradient and consequently 
the /^-gradient. This hypothesis leads to a diffusion coef- 
ficient D sll car given by 



D, 



shear 



= 4 



K 



1 / Jhfi 

-a 0.8836 ft— 

4 \ d In r 



-(V'-V) 



.(16) 



The term V' — V in Eq. (16) expresses either the stabilizing 
effect of the thermal gradients in radiative zones or its 
destabilizing effect in semiconvective zones (if any). When 
the shear is negligible, Z3 s hcar tends towards the diffusion 
coefficient for semiconvection by Langer et al. ( |1983 ) in 
semiconvective zones. When the thermal losses are large 
(V' = V), it tends towards the value 



a{K/Nl A ) 0.8836 fl 



dlnr 



(17) 



given by Zahn (1992). Eq. (16) is completed by the three 



following equations expressing the thermal effects (Maeder 
1997] ) 



-Dshcar = 2.KT 



V - V = 



V = 



rad 



r + i 



(v 



ad 



V). 



(18) 



(19) 



The system of 4 equations given by Eqs. (16), (18) and 
(19) form a coupled system with 4 unknown quantities D, 
r, V and V. The system is of the third degree in T. When 
it is solved numerically, we find that as a matter of fact 
the thermal losses in the shears are rather large in massive 
stars and thus that the Peclet number Pe is very small 
(of the order of 10~ 3 to 10~ 4 , see Sect. 4.2). For very low 
Peclet number Pe = 6T, the differences (V — V) are also 
very small as shown by Eq. (19). Thus, we conclude that 
Eq. (16) is essentially equivalent, at least in massive stars, 
to the original Eq. (17) above, as given by Zahn (1992). 
We may suspect that this is not necessarily true in low and 
intermediate mass stars since there the Pe number may be 
larger. Of course, the Reynolds condition -D s hcar > ^vRe c 
must be satisfied in order that the medium is turbulent. 
The quantity v is the total viscosity (radiative + mole- 
cular) and Re c the critical Reynolds numb er esti mate d 
to be around 10 (cf. Denissenkov et al. 1999 Zahn 1992 ). 
The numerical results in Sect. 4 will show the values of the 
various parameters and also indicate that the conditions 
for the occurence of turbulence are satisfied. 

The physical treatment around the core also depends 
on the choice of the criterion for convection. In current 
literature, there are at least three basic sets of assump- 
tions: a) Ledoux criterion, which leads to small cores, 
b) Schwarzschild's criterion, which gives "medium size" 
cores, c) Schwarzschild's criterion and overshooting, which 
gives large cores. Here, we choose to apply the intermedi- 
ate solution b), which means that the above equation 16 is 
only applied in fully radiative regions. Some consequences 
of this choice of convective criterion are discussed in Sect. 
3 below. 

2.9. Initial compositions, opacities, nuclear reactions and 
other model ingredients 

For purpose of comparison, we adopted here the same 
physical ingredients as for the solar metallicity models 
computed by Meynet et al. (1994). The only exceptions, 
apart from the inclusion of the effects induced by ro- 
tation described above, are the following: we use the 
Schwarzschild criterion for convection without overshoot- 
ing and the mass loss rates are as indicated in Sect. 2.4. 

3. The non rotating stellar models 

For comparison purpose, we have computed non-rotating 
stellar models with exactly the same physical ingredi- 
ents as the rotating ones. The corresponding evolutionary 
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tracks and lifetimes are presented in Fig. 8 and Table 1. 
These stellar models have similar properties as older ones 
computed with the Schwarzschild criterion for convection. 
For instance, our non-rotating 9 M Q stellar model has 
similar H- and He-burning lifetimes as the 9 M Q model 
computed with the Schwarzschild criterion by Bertelli et 
al. (1985). Our models also well agree with more recent 
computations. Indeed our MS lifetimes are similar within 
about two percents to the ones obtained by Heger et al. 
(2000) for their non-rotating stellar models. These authors 
used the Ledoux criterion with semiconvection. However, 
since during the MS phase the convective core mass de- 
creases, one expects, for this phase, only small differences 
between the models computed with the Ledoux and the 
Schwarzschild criterion. These comparisons show that the 
present non-rotating stellar models well agree with results 
obtained with different stellar evolutionary codes. 

Many papers (e.g. Bertelli et al. 1985; Maeder & 
Meynet 1989; Chin & Stothers 1990; Langer & Maeder 
1995; Canuto 2000) have discussed the effects of different 
criteria for convection on massive star evolution. Since our 
main aim in this paper is to emphasize the effects of ro- 
tation, we shall restrict ourselves to briefly mention some 
differences with previous grids of models computed by the 
Geneva group. 

In the present non-rotating stellar models, computed 
without overshooting, the convective core masses are 
smaller than in the models by Schaller et al. (1992) and 
Meynet et al. (1994). As a consequence, the stars evolve 
at smaller luminosities. The turn-off point at the end of 
the MS occurs for higher T e ff and lower luminosities, re- 
ducing the extent of the MS width. The MS lifetimes are 
reduced. For initial masses below about 25 M Q , the He- 
burning lifetime is increased implying an increase of the 
ratio tjjo/iH of the lifetimes in the He- and H-burning 
phases. For the higher initial masses, the effects of the 
stellar winds become very important and dominate the 
effects due to a change of the criterion for convection. 

As a numerical example, in our non-rotating 20 M 
stellar model, the convective core mass is reduced by 1— 
1.2 M Q during the whole H-burning phase compared to 
its value given by Meynet et al. (1994). At the end of the 
MS phase, the mass of the helium core has decreased by 
~ 20% with respect to its value in models with a mod- 
erate overshoot. The position of the turn-off point has a 
logL/L Q decreased by 0.05 dex and a logT ff increased by 
0.04 dex. The MS lifetime is decreased by about 10%, the 
ratio tji e /tn passes from 10% in the models of Meynet et 
al. (1994) to 14% in the present case. 

4. The evolution of the internal rotation law f2(r) 

4-1- The initial convergence of the rotation law 

Let us first mention that even in some recent works the as- 
sumption of solid body rotation is often considered. How- 
ever, it is much more physical to examine the evolution of 
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Fig. 1. Initial evolution of the angular velocity fl as a 
function of the distance to the center in a 20 M Q star 
with v m i = 300 km s _1 . Starting from a flat profile, the 
solutions rapidly converge towards an equilibrium solu- 
tion. The time interval between two consecutive curves is 
19 200 years. 

H(r) resulting from the transport of angular momentum 
by shears, meridional circulation and convection, from the 
effects of central contraction and envelope expansion and 
from the losses of angular momentum by stellar winds. 
In particular, the large losses of angular momentum at 
the surface lead to a redistribution of the angular mo- 
mentum in the interior by meridional circulation, shear 
turbulent diffusion and convection. The whole problem is 
self-consistent, since the meridional circulation and shear 
transport depend in turn on the degree of differential ro- 
tation. Thus, the full system of equations has to be solved 
in order to provide the internal evolution of fi. 

Fig. H shows the results for the initial convergence on 
the zero age main sequence of f2(r) in a 20 M Q model 
with solar composition. The initial model on the zero-age 
main sequence is supposed to rotate uniformly and with 
an initial surface velocity of 300 km s _1 . Very short 
time steps of the order of 1920 yr are taken in the initial 
calculations for obtaining the internal equilibrium rota- 
tion. The results are shown in Fig. [j] for every 10 th model, 
i.e. with time intervals of 19 200 yr. We see the very fast 
initial changes of O(r), mainly due to the meridional cir- 
culation. The circulation velocity, which is very large, is 
also positive everywhere at the beginning (cf. Fig. ||). This 
means that the circulation rises along the polar axis and 
goes down at the equator, thus bringing angular momen- 
tum inwards. As a consequence the angular velocity of the 
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Fig. 2. Same as Fig. [j] for U(r) the vertical component 
of the meridional circulation velocity (cf. Eq. 4). Time 
proceeds from top to bottom. 

core increases, while that in the envelope goes down. It is 
important to realize that here these fast changes of the 
rotation law are not due, as later in the evolution, to core 
contraction and envelope expansion. 

We see that the changes are big at the beginning and 
then smaller. The profiles of Q(r) rapidly converges to- 
wards an equilibrium-profile in a time of about 1 to 2% of 
the MS lifetime tyis ■ This is in agreement with t he re sults 
for 10 and 30 M stars by Denissenkov et al. ( |l999|) . As 
noted by these authors, the timescale for the adjustment of 
the rotation law, i.e. t c ; rc ~ R/U , is short with respect to 
the MS lifetime for any significant rotation. Later on dur- 
ing the bulk of MS evolution, the profiles Q(r) will evolve 
more slowly. As emphasized in Sect. 2.6, this timescale i c ; rc 
is not characteristic for the mixing of chemical elements. 
We also notice that the degree of differential rotation in 
the equilibrium profile is rather modest (cf. Zahn, 1992). 

Fig. |^ shows the corresponding initial evolution of the 
vertical component U(r) of the velocity of the meridional 
circulation in the 20 M Q model. As seen above, initially 
this velocity is large and positive everywhere, which ex- 
plains the fast evolution of f2(r) in Fig. [l]. Then, U(r) de- 
creases and becomes negative in the very external layers. 
The physical reason for that is the so called Gratton-Opik 
term of the form — 2 - K Gp wmc h i s contained in the expres- 
sion of Eq in Eq. (12). When the density is very low, as in 
the outer regions, this negative term becomes important 
and produces an inverse circulation. This means that the 
circulation has two big cells, an internal cell rising along 
the polar axis and an external one descending at the pole. 



Fig. 3. Non-stationary (continuous line) and stationary 
(broken line) solutions for U(r) the vertical component of 
the velocity of the meridional circulation (cf. Eq. 4) as a 
function of the distance to the center in the 20 M model 
at the end of the initial convergence. 

The velocity U (r) also converges towards an equilibrium 
distribution, characterized by small veloci ties U (r), a re- 
sult in agreement with Denissenkov et al. (199E). 

It is interesting to compare the stationary solution of 
U{r) given by Eq. (5) and the non-stationary solution 
given by Eq. (12). This is done in Fig. |^, where the curve 
in continuous line corresponds to the asymptotic distribu- 
tion reached in the non-stationary regime illustrated in 
Fig. ||. This curve shows positive values of U(r) in the in- 
ner regions and negative values in the outer ones. On the 
contrary, the solution (broken line) obtained in the sta- 
tionary approximation given by Eq. (5) is always positive. 
This is a logical consequence of the approximation made: 
as f2(r) decreases outwards, only positive values of U(r) 
are possible. Said in other words, the outward transport 
by shears can only be compensated by the inward trans- 
port by meridional circulation, i.e. by positive values of 
U{r). 

However, in non-stationary models, this compensation 
is not achieved since the two curves in Fig. [| are very dif- 
ferent. One concludes that the stationary approximation 



which has been used in some works (cf. Urpin et al. 1996) 
is not satisfactory. 

The stationary solution can also be viewed as giv- 
ing a velocity which is the inverse of a velocity which 
could be associated to the diffusive transport by shears. 
From the comparison of the two curves in Fig. |^, one thus 
immediately concludes that the meridional circulation is 
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much more efficient for the transport of the angular mo- 
mentum than the shear instabilities. The justification for 
this claim is the following one. In the stationary situa- 
tion, shear transport and meridional circulation compen- 
sate each other and are thus of the same magnitude. The 
fact that U (r) in the full non-stationary case treatment is 
much larger than in the stationary case as shown in Fig. 3 
implies that the effects of meridional circulation on the 
transport of angular momentum are much larger than the 
effects of the shear. This conclusion can also be deduced 
from the comparison of the timescales for the two pro- 
cesses. In general one has that i c irc = § < £ 8 hear = jj^ — ■ 
For the chemical elements, the transport by the shear 
instabilities is more efficient than the transport by the 
meridional circulation as discussed in Sect. 2.6 where we 
have seen that Z? ff is reduced by horizontal turbulence. 

4-2. Internal evolution of Cl(r) 

The evolution of f2(r) during the MS evolution of a 20 
Mq star is shown in Fig. |I| We notice that initially (i.e. 
for high values of the central H-content X c ), the degree 
of differential rotation is small, then, it grows during the 
evolution. The ratio between the central and surface value 
of f2 remains however small during the MS evolution. The 
rotation rate does not vary by more than about a factor 
two throughout the star. Although this degree of differen- 
tial rotation is not large, it plays an essential role in the 
shear effects and the related transports of angular momen- 
tum and of chemical elements. The same remark also ap- 
plies to the effects of differential rotation on the transport 
by meridional circulation, since U(r) critically depends on 
the derivatives of VL in the star (cf. Zahn 1992; Maeder & 
Zahn 1998). 

We notice a general decrease of £l(r), even in the con- 
vective core which is contracting. The main reason is mass 
loss at the stellar surface, which removes a substantial 
fraction of the total angular momentum. This makes Q(r) 
to decrease with time everywhere, because of the inter- 
nal transport mechanisms, which ensure some coupling 
of rotation. The same behaviour was obtained by Langer 
(1998) in the case of rigid rotation (i.e. in the case of max- 
imum coupling). In the present models, shear transports 
the angular momentum outwards, which reduces the core 
rotation, while circulation makes an inward transport in 
the deep interior and an outwards transport in the exter- 
nal region. This is responsible for the progressive flatten- 
ing of the curves above r = 4i? Q in Fig. |]. 

From the end of the MS evolution onwards, i.e. when 
X c < 0.05, central contraction becomes faster and starts 
dominating the evolution of the angular momentum in 
the center. The central f2 grows quickly (cf. dotted curve 
in Fig. and this will in principle continue in the fur- 
ther evolutionary phases until core collapse. The evolu- 
tion of the angular momentum after the H-burning phase 
(X c = 0) is mainly dominated by the local conservation of 




Fig. 4. Evolution of the angular velocity as a function 
of the distance to the center in a 20 M star with Vi n i 
= 300 km s _1 . X c is the hydrogen mass fraction at the 
center. The dotted line shows the profile when the He- 
core contracts at the end of the H-burning phase. 



the angular momentum, since the secular transport mech- 
anisms have little time to operate. This is the assumption 
we are making in the present models. During these phases, 
the chemical elements continue to be rotationally mixed 
in the radiative layers mainly through the effect of the 
shear. In the further evolutionary phases, the strong cen- 
tral contraction will lead to very large central rotation, 
unless some other processes as fast dynamical instabilities 
or magnetic braking occur in these stages. 

The evolution of U(r) during MS evolution is shown 
m Fig. |. The outer zone with inverse circulation progres- 
sively deepens in radius during MS evolution due to the 
Gratton-Opik term, because a growing part of the outer 
layers has lower densities. Also, we notice that the velo- 



cities U(r) are small in general (cf. Urpin et al. 1996) and 
tests have shown us that, contrary to the classical result of 
the Eddington-Sweet circulation (cf. Mestel 1965), U(r) 
depends rather little on the initial rotation. 

The deepening of the inverse circulation has the conse- 
quence that the stationary and non-stationary solutions 
differ more and more as the evolution proceeds, since as 
said above no inverse circulation is predicted by the sta- 
tionary solution. Thus we conclude that the stationary 
solutions are much too simplified. Also, in Fig. we see 
that the values of U (r) become more negative in the outer 
layers for the model at the end of the MS phase, when 
central contraction occurs. Again, this effect is due to 
the Gratton-Opik term. The results of Fig. 5 also shows 
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Fig. 5. Same as Fig. [| for U(r). 

that the application of a diffusive treatment to meridional 
circulation transport is unappropriate. Paradoxically, the 
problem is less serious in regions which exhibit an inverse 
circulation, since there the signs of the effects are at least 
the same in the two treatments (however even there the 
equations for advection and diffusion arc different). 

The various diffusion coefficients inside a 20 M Q star 
are shown in Fig. 6, when the hydrogen mass fraction X c 
at the center is equal to 0.20. We notice that in general 
K > Dh > -Dshcar > D c s > v. Following the thermal 
diffusivity K, the coefficient of horizontal diffusion Dh is 
the second largest one, this is consistent and necessary for 
the validity of the assumption of shellular rotation (Zahn, 
1992|) Since we have Ashear = 2KT according to Eq. (18), 
Fig. |6| shows that T is typically of the order of 10~ 3 to 
10 -4 . The coefficient D e g, expressing the transport of the 
chemical elements by the meridional circulation with ac- 
count of the reduction by the horizontal turbulence, is only 
slightly smaller (by a factor of about 3) than -D s hcar in the 
interior. One must be careful that D e g concerns only the 
tranport of the elements and not that of angular momen- 
tum, for which meridional circulation is the largest effect. 
Fig. 6 also shows the values of the total viscosity is, and 
we notice that consistently the Prandtl number v/ K is of 
the order of 10~ 5 to 10~ 6 . 



5. HR diagram, lifetimes and age estimates 

5.1. The Main-Sequence evolution 

Evolutionary tracks of 20 M© models at solar metallicity 
for different initial velocities are plotted on Fig. @. The 
T e ff for a rotating star has been defined in Sect. 2.3. On 



Fig. 6. Internal values of K the thermal diffusivity, Dh 
the coefficient of horizontal turbulence, D s hear the shear 
diffusion coefficient, D c g the effective diffusivity (cf. Eq. 9) 
and v the total viscosity (radiative + molecular) in the 
radiative envelope of a 20 Mq star with an initial v- lv \ = 
300 km s _1 . The lagrangian mass coordinate is given on 
the upper scale. Here, the hydrogen mass fraction at the 
center X r = 0.20. 



and near the ZAMS, rotation produces a small shift of 
the tracks towards lower luminosities and T g. This effect 
is due to both atmospheric distorsions and to the low- 
ering of the effective gravity (see e.g. Kippenhahn and 
Thomas 1970 ; Maeder and Peytremann 1970 Collins and 
Sonncborn 1977 ). At this stage the star is still nearly ho- 
mogeneous. When evolution proceeds, the tracks with ro- 
tation become more luminous than for no rotation. This 
results from essentially two effects. On one side, rotational 
mixing brings fresh H fuel into the convective core, slow- 
ing down its decrease in mass during the MS. For a given 
value of the central H-content, the mass of the convective 
core in the rotating model is therefore larger than in the 
non-rotating one and thus the stellar luminosity is higher 
(Maeder |1987[ Talon et al. ~ 



19971 Heger et al 
= 300 km s 



200q ). As 
1 models 



a numerical example, in the ■ 
the He-cores at the end of the MS are about 20% more 
massive than in their non-rotating counterparts. This is 
equal to the increase obtained by a moderate overshooting 
(see Sect. 3). On the other side, rotational mixing trans- 
ports helium and other H-burning products (essentially 
nitrogen) into the radiative envelope. The He-enrichment 
lowers the opacity. This contributes to the enhancement 
of the stellar luminosity and favours a blueward track. In- 
deed, in Fig. 0, one sees that when the mean velocity on 
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velocity is v- u 



300 km s . There is little difference 
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Fig. 7. Evolutionary tracks for rotating 20 M© models 
with different initial velocities. The mean equatorial ve- 
locities v during the MS are indicated. The corresponding 
initial velocities can be seen in Table 2. 

the MS becomes larger than about 130 km s _1 , the T e g 
at the end of the MS increases when rotation increases. 

For sufficiently low velocities, rotation acts as a small 
overshoot, extending the MS tracks towards lower T e g. 
This results from the fact that, at sufficiently low rota- 
tion, the effect of rotation on the convective core mass 
overcomes the effect of helium diffusion in the envelope. 
Indeed, for small rotation, the time required for helium 
mixing in the whole radiative envelope is very long, while 
hydrogen just needs to diffuse through a small amount of 
mass to reach the convective core. 

Since not all the stars have the same initial rotational 
velocity, one expects a dispersion of the luminosities at the 
end of the MS. For the 20 M© models shown on Fig. [?] one 
sees a difference of ~ 0.3 in Mboi between the luminosities 
of the low and fast rotators at the end of the MS. Rota- 
tion induces also a scatter of the effective temperatures at 
the end of the MS. In reality, the dispersion results from 
both different initial velocities and also, for a given initial 
velocity, from different angles between the axis of rotation 
and the line of sight. Indeed due to the von Zeipel theorem 
( |1924D the star appears bluer seen pole-on than equator- 
on. When integrated over the visible part of the star, the 
effects due to orientation can reach a few tenths of a mag- 
nitude in luminosity and a few hundredths in logT c ff (cf. 
Maeder and Peytreman 1970). 

Figure ^| shows the evolutionary tracks of non-rotating 
and rotating stellar models for initial masses between 9 
and 120 M©. For the rotating stellar models, the initial 



between tracks with 
Talon et al 



1997) 



= 200 or 300 km s" 1 (see also 



If the effects behaved like f ini , there 
would be larger differences. The present saturation effect 
occurs because outward transport of angular momentum 
by shears are larger when rotation is larger, also a larger 
rotation produces more mass loss, which makes a larger 
reduction of rotation during the evolution. Let us note 
however that for some surface abundance ratios as N/C 
or N/O (see Table 1), the increase from Vi n [ = 200 to 300 
km s _1 produces significant changes. Thus, the similarity 
of the evolutionary tracks does not necessarily imply the 
similarity of the surface abundances for these elements. 

Rotation reduces the MS width in the high mass range 
(M <; 40 M ). Let us recall that when the mass increases, 
the ratio of the diffusion timescale for the chemical ele- 



ments to the MS lifetime decreases (Maeder 1998). As a 
consequence, starting with the same Ui n i on the ZAMS, 
massive stars will be more mixed than low mass stars at 
an identical stage of their evolution. This reduces the MS 
width since greater chemical homogeneity makes the star 
bluer. Moreover, due to both rotational mixing and mass 
loss, their surface will be rapidly enriched in H-burning 
products. These stars will therefore enter the Wolf-Rayet 
phase while they are still burning their hydrogen in their 
core. This again reduces the MS width. For initial masses 
between 9 and 25 M©, the MS shape is not much changed 
by rotation at least for Vi n [ < 300 km s _1 . 

5.2. The post-Main-Sequence evolution 

The post-MS evolution of the most massive stars (M > 40 
Mq) which become W-R stars will be discussed in a forth- 
coming paper. We shall just mention one point of general 
interest here: for low or moderate rotation, the convec- 
tive core shrinks as usual during MS evolution, while for 
high masses (M <; 40 M ) and large initial rotations 
(fp— > 0.5), the convective core grows in mass during 
evolution. This latter situation occurs in the fast rota- 
ting 60 Mq model shown on Fig. || These behaviours, i.e. 
reduction or growth of the core, determine whether the 
star will follow respectively the usual redwards MS tracks 
in the HR diagr am, o r whethe r it w ill bifurcate to the 
blue (cf. Maeder 1987 ; Langer 1992 ) towards the classi- 
cal tracks of homogeneous evolution (Schwarzschild 1958 ) 
and likely produce W-R stars. 

The stars with initial masses between 15 and 25 M© 
become red supergiants (RSG). Rotation does not change 
qualitatively this behaviour but accelerates the redwards 
evolution, especially for the 15 and 20 M Q models. As a 
numerical example, for an initial i>i n j = 300 km s _1 , the 
model stars burn all their helium as red supergiants at 
T e g below 4000 K, while the non-rotating models spend a 
significant part of the He-burning phase in the blue part 
of the HR diagram: for the non-rotating 15 and 20 M© 
models, respectively 25 and 20% of the total He-burning 
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Fig. 8. Evolutionary tracks for non-rotating (dotted lines) and rotating (continuous lines) models with solar metallicity. 
The rotating models have an initial velocity i>i n i of 300 km s _1 . For purpose of clarity, only the first part of the tracks 
for the most massive stars (M > 40 M©) is shown. Portions of the evolution during the W-R phase for the rotating 
massive stars are indicated by short-dashed lines. The long-dashed track for the 60 M model corresponds to a 
very fast rotating star (i; ini ~ 400 km s _1 ), which follows a nearly homogeneous evolution. Only the beginning of its 
evolution is shown. 



lifetime is spent at logT c ff > 4.0. The behaviour of the 
rotating models results mainly from the enhancement of 
the mass loss rates. This effect prevents the formation of 
a big intermediate convective zone and therefore favours a 
rapid evolution toward the RSG phase (Stothcrs and Chin 



1979; Maeder 1981). Let us note that the dispersion of the 



initial rotational velocities produces a mixing of the above 
behaviours. 

Very interestingly for the 12 M Q model a blue loop 
appears when rotation is included. This results from the 
higher luminosity of the rotating model. The higher lu- 
minosity implies that the outer envelope is more ex- 
tended, and is thus characterized by lower temperatures 
and higher opacities at a given mass coordinate. As a con- 
sequence, in the rotating model during the first dredge-up, 
the outer convective zone proceeds much more deeply in 



mass than in the non-rotating star. Typically in the non- 
rotating model the minimum mass coordinate reached by 
the outer convective zone is 6.6 M© while in the rotating 
model it is 2.6 M . This prevents temporarily the exten- 
sion in mass of the He-core and enables the apparition of a 
blue loop. Indeed the lower the mass of the He-core is, the 
lower its gravitational potential. According to Lauterborn 



et al. (1971, see also the discussion in Maeder and Meynet 



1989), a blue loop appears when the gravitational poten- 



tial of the core $ c is inferior to a critical potential $crit 
depending only on the actual mass of the star which is 
about the same for the rotating and non-rotating model. 
This explains the appearance of a blue loop in the 12 M© 
rotating model. For the 9 M© model, the minimum mass 
coordinate reached by the outer convective zone is not 
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Table 1. Properties of the stellar models at the end of the H-burning phase, at the blue supergiant (BSG) or LBV 
stage (see text) and at the end of the He-burning phase. The masses are in solar mass, the velocities in km s _1 , the 
lifetimes in million years and the abundances in mass fraction. 
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much affected by rotation and the models with and with- 
out rotation present very similar blue loops. 



5.3. Masses and mass-luminosity relations 

When rotation increases, the actual masses at the end of 
both the MS and the He-burning phases become smaller 
(cf. Tables |l| and |J). Typically the quantity of mass lost 
by stellar winds during the MS is enhanced by 60-100% in 
rotating models with v lnl — 200 and 300 km s _1 respec- 
tively. For stars which do not go through a Wolf-Rayet 
phase, the increase is due mainly to the direct effect of 
rotation on the mass loss rates (in the present models 
through the formula proposed by Friend and Abbott 1986 ) 
and to the higher luminosities reached by the tracks com- 
puted with rotation. The fact that rotation increases the 
lifetimes also contributes to produce smaller final masses. 
For the most massive stars (M > 60 M ), the present ro- 
tating models enter the Wolf-Rayet phase already during 
the H-burning phase (see also Maeder 1987; Fliegner & 
Langer 1995; Meynet 1999, 2000b). This reduces signifi- 
cantly the mass at the end of the H-burning phase. 

As indicated in Sect. 5.1, the initial distribution of the 
rotational velocities implies a dispersion of the luminosi- 



ties at the end of the MS. This effect introduces a sig- 
nificant scatter in the mass-luminosity relation (Langer 
1992; Meynet 1998 ), in the sense that fast rotators are 
overluminous with respect to their actual masses. This is 
especially true in the high mass star range in which the 
luminosity versus mass relation flattens. This may explain 
some of the discrepancies between the evolutionary masses 
and the direct mass estimates in some binaries (Penny et 
al. p999h. 



Let us end this section by saying a few words about the 
mass discrepancy problem (see e.g. Herrero et al. 2000). 
For some stars, the evolutionary masses {i.e. determined 
from the theoretical evolutionary tracks) are greater that 
the spectroscopically determined masses. Interestingly, ac- 



cording to Herrero et al. (2000), only the low gravity ob- 
jects present (if any) a mass discrepancy. Even if most of 
the problem has collapsed and was shown to be a result of 
the proximity of O-stars to the Eddington limit (Lamers 
and Leitherer 1993; Herrero et al. 1999) and of the large 
effect of metal line blanketing not usually accounted for 
in the atmosphere models of massive stars (Lanz et al. 
1996), some discrepancy seems to be still present. The re- 



maining mass discrepancy may arise, in part, from the use 
of non-rotating models for determining the evolutionary 
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masses. High rotation produces larger He-cores and He- 
rich envelopes, all that implies overluminous stars with 
lower gravity as evolution proceeds. This can occur even 
for the slow rotating objects because they could have had 
a sufficiently high initial rotational velocity. The observa- 
tion that the mass discrepancies are found only for low 
gravity objects may reflect the fact that rotation implies 
more and more important changes in the logoff versus 
log Teg plane when evolution proceeds. Indeed, looking at 
Fig. |l6| below one sees that at high gravity (i.e. at an 
early evolutionary stage), one does not expect any mass 
difference when using rotating or non-rotating tracks. In 
contrast, the mass difference between the rotating and the 
non-rotating tracks become more and more important in 
the high mass star range and for the low values of logoff- 
The rotating 40 M Q track crosses the non-rotating 60 M Q 
model on Fig. [l6| indicating that differences of ~30% in 
mass are quite possible. 

For the fast rotating objects, a part of the mass dis- 
crepancy may also be due to a possible underestimate of 
the gravity (see Sect. 7.1). Indeed, an underestimate of the 
gravity would also imply an underestimate of the spectro- 
scopically determined masses (Herrero et al. 2000). 

5.4- Lifetimes and isochrones 

Table |l| presents some properties of the models. Column 1 
and 2 give the initial mass and the initial velocity re- 
spectively. The mean equatorial rotational velocity v dur- 
ing the MS phase is indicated in column 3. This quantity 
is defined by 

v = l/t H / v(t)dt, 
Jo 

where tu is the duration of the H-burning phase given 
in column 4, except for those stars which enter the Wolf- 
Rayet phase while still burning their hydrogen in their 
core, i.e. for the rotating 60 and 120 M models, for which 
t H has been replaced by the duration of the O-type star 
phase. The H-burning lifetimes the masses M, the 
equatorial velocities v, the helium surface abundance Y s 
and the surface ratios (in mass) N/C and N/O at the end 
of the H-burning phase are given in columns 4 to 9. The 
columns 10 to 13 present some properties of the models 
when the star is a blue supergiant (BSG) or an LBV star. 
For stellar models with M< 40 M , the "BSG stage" in 
Table [j] corresponds either to the stage when log T c $ — 
4.0 during the first crossing of the Hertzsprung-Russel 
diagram or to the bluest point on the blue loop if any 
for M< 12 M Q . For non-rotating stellar models with M 
> 60 M , the "LBV stage" corresponds to the point when 
the star has lost half of the matter ejected by the stellar 
winds between the end of the MS and the entrance into 
the W-R phase. In the case of rotating models, the "LBV 
stage" corresponds to the period during which the surface 
velocity becomes critical and huge mass loss rates ensue. 



Table 2. Properties of 20 M Q models at the end of the 
MS for different initial velocities. The velocities are in km 
s _1 , the lifetimes in million years, the masses in solar mass 
and the abundances in mass fraction. 
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Again, here we choose a model in the middle of this phase. 
The columns 14 to 19 present some characteristics of the 
stellar models at the end of the He-burning phase and tu e 
is the He-burning lifetime. In Table || some properties of 
20 M Q models at the end of the MS are indicated, v- ln i and 
v have the same meaning as above. 

From Table one sees that for Z = 0.020 the lifetimes 
are increased by about 20-30% when the mean rotational 
velocity on the MS increases from to -200 km s" 1 . This 
modest increase is explained by the fact that even if there 
is more fuel available in the core, the luminosity is also 
increased. From the data presented in Table ^, one can 
deduce a nearly linear relation between the relative en- 
hancement of the MS lifetime, A% and v, where 

AiH(U) = Mu)-< h (0)]Ah(0). 

One obtains, 



where v is in km s _1 . This relation reproduces the values 
of Atff(u)/tjj(0) from table || with an accuracy better 
than 5%. It also applies with the same accuracy to the 
values listed in Table |l| for the masses between 15 and 40 
Mq . The He-burning lifetimes are less affected by rotation 
than the MS lifetimes. The changes are less than 10%. The 
ratios tne/tH of the He to H-burning lifetimes are only 
slightly decreased by rotation and remain around 10-15%. 

On Fig. ^| isochrones for ages between about 8 and 20 
10 6 yr, computed from non-rotating and rotating stellar 
models are presented. The "rotating" isochrones are com- 
puted from the models with an initial rotational velocity 
v ini of 200 km s" 1 on the ZAMS. 

At a given age, the upper part of the "rotating" 
isochrones are bluer and more luminous. Typically, for an 
age equal to about 20 • 10 6 years (log age =7.3) the red- 
dest point on the MS is shifted by 0.015 dex in logT c ff and 
by - 0.6 in Mboi when rotation is taken into account. An 
isochrone with rotation is almost identical to an isochrone 
without rotation with log age smaller by 0.1 dex. This has 
for consequence that rotation slightly increases the age 
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Fig. 9. Isochrones computed from stellar evolutionary 
tracks for the solar metallicity. The dashed and continuous 
lines correspond to the case of non-rotating and rotating 
stellar models respectively. In this last case, models have 
an initial velocity of 200 km s _1 . The logarithms of 
the ages (in years) label the isochrones computed from the 
models with rotation. 

associated to a given cluster. From Fig. 9, one sees that 
the "rotating" isochrone for an age equal to 20 • 10 6 years 
has the same luminosity at the turn-off than the "non- 
rotating" isochrone for an age equal to 16- 10 6 years. Thus 
rotation increases the age estimate by about 25%. We may 
wonder whether this effect explains the age difference be- 
tween the estimates based on the upper MS and the esti- 
mates based on the lithium content of the very low mass 
stars (see e.g. Martin et al. 1998; Barrado y Navascues et 
al. 1999). One must also account for the dispersion of ro- 
tational velocities and possibly of the orientation angles. 
These two effects introduce some dispersion in the way 
stars are distributed in the HR diagram and thus affect 
the interpretation of the clusters' observed sequences (cf. 
Maeder |1971| ). 

If a bluewards track occurs, as for very massive stars 
with fast rotation, the larger core and mixing lead to much 
longer lifetimes in the H-burning phase. In this case, the 
fitting of time-lines becomes hazardous. 

6. Evolution of the rotational velocities 

6.1. Model results for stars with a large mass loss 

Figs. 10 and |ll| show the evolution of v and of j^— as 
a function of the age for the present models. Firstly, we 
notice that for models without mass loss, as shown for the 



Fig. 10. Evolution of the surface equatorial velocity as a 
function of time for stars of different initial masses with 
t>ini = 300 km s _1 . The continuous lines refer to solar 
metallicity models, the dotted line corresponds to a 20 
M star with Z = 0.004. The dashed line corresponds to 
a 20 M star without mass loss. 



20 M with M = 0, v and go up fastly so that the 
critical velocity would be reached near the end of the MS 
phase. The current model of 20 M with mass loss show 
a significant decrease of v, while the critical ratio remains 
almost constant during most of the MS phase. Figs. [r] 
and [ll] show how fastly rotation decreases at the surface 
of the most massive stars, which lose a lot of mass. Con- 
sistently we see that the reduction of the surface rotation 
is much larger for the more massive stars. This is of course 
a consequence of the removal of large amounts of angular 
momentum by the stellar winds. The effect is amplified by 
the increase of the mass loss in fast rotators (Eq. 2). We 
see that the decrease of jp-^ is so strong that it will pre- 
vent a massive star to reach the critical velocity near the 
end of the MS phase. If the star makes extended excur- 
sions in the HR diagram at the end of the MS like is the 
case for the 60 M model, then it may reach the critical 
velocity. The specific case of stars close to the fi-limit will 
be examined in a future study, since this requires some 
further theoretical developments. 

Let us mention here that the present results differ 



from those obtained by Sackmann and Anand (197C) and 
Langer (1997, 1998] ). Indeed, these authors find that the 
star reaches the break-up limit during the MS phase. As 
an example, in a 60 M star, even a model with an ini- 
tial Vi n i of 100 km s _1 reaches the break-up limit near 
the end of the MS-phase (Langer 1998). This result leads 
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Fig. 11. Same as Fig. |l^ for the ratio f2/f2 C rit of the angu- 
lar velocity to the break-up velocity at the stellar surface. 



effects of rotationnally induced mixing. We see here that 
modifying these hypothesis (and also using other prescrip- 
tions for the mass loss rates) lead to very different results. 
One of the first step to clarify the situation is to deter- 
mine which expression for u cr it (cf. Glatzel 1998; Langer 
1998) is the correct one and how rotation affects the mass 
loss rates. These developments, now in progress, will be 
particularly needed for the study of the evolution of the 
most massive stars, like a 120 M Q model, which have a 
high value of the Eddington factor. Also this is important 
for the formation and evolution of W-R stars, which will 
be studied in a further work. 

Fig. |l2] shows the evolution of v with age for models of 
a 20 M star with different initial velocities w; n i from 50 
to 580 km s _1 . We see that the decrease in the surface v is 
larger for larger initial rotation. This is a consequence of 
the larger mass loss rates in fast rotators. We notice some 
convergence (cf. Langer 1998) of the curves for the large 
initial uj^. This convergence would be more pronounced 
and affect the models with lower if the dependence 
of M vs v would be stronger than given by Eq. (2). This 
certainly reduces the scatter of v near the end of the MS 
phase, but does not produce a full convergence. 
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Fig. 12. Evolution of the surface equatorial velocity as 
a function of time for 20 Mq stars with different initial 
velocities. 

Langer to conclude that most massive stars may reach 
the break-up limit or the so-called fi— limit during their 
MS evolution. Let us however emphasize here that such a 
conclusion is based on a particular definition of u cr it still 
subject to discussion (see Sect. 2.4), on the assumption of 
solid body rotation, and on models not accounting for the 



6.2. Results for stars with lower mass loss 

The models of 12 and 15 M Q show curves in Figs. [ic| 
and [ll] with little reduction of v, while there are slight 
increases of the critical ratio jp— during MS evolution. 
A peak is reached during the overall contraction phase at 
the end of the MS phase. Then, v and the critical ratio 
go down as the star moves to the red supergiant phase. 
During such a fast evolution, we may say that the rotation 
evolves almost like the case of simplified models, where the 
angular momentum is conserved locally. The large growth 
of the radius just implies a decrease of v and of . Later 
during the blue loops, where the Cepheid instability strip 
is crossed, the rotation velocity becomes very large again 
and could easily become close to critical. This behaviour, 
also found by Heger & Langer (1998), results from the 
stellar contraction which concentrates a large fraction of 
the angular momentum of the star (previously contained 
in the extended convective envelope of the RSG) in the 
outer few hundredths of a solar mass. This result suggests 
that rotation may also somehow influence the Cepheid 
properties, in addition to the increase of stellar luminosity 
discussed above. 

Amazingly, we notice that the stars with initially low 
mass loss during the MS phase, like for stars with M < 12 
Mq, have more chance to reach the break-up velocities 
and thus huge mass loss than the more massive stars which 
lose a lot of mass on the MS. It is somehow surprising that 
little mass loss during the MS may favour large mass loss 
rates at the end of the MS phase. This may explain why 
stars close to break-up, like the Be stars, do not form 
among the O-type stars, but mainly among the B-type 
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Fig. 13. Evolution of the equatorial surface velocities along the evolutionary tracks in the HR diagram starting from 
Wini = 300 km s _1 . For purpose of clarity, only the first part of the 40 M Q track is shown. 



stars, where we see that the ratio may increase during 
the MS phase. Another related observation is the fact that 
the relative number of Be-stars with respect to B-type 
stars is much higher in t he LM C and SMC than in the 
Galaxy (cf. Maeder et al. |l999|) . In the LMC and SMC, 
due to the lower metallicity, the average mass loss rates 
are lower and thus these stars may keep higher rotation 
in general and thus form more Be stars. This explanation 
does not exclude differences in the distribution of v- m i as 
well. 



6.3. The rotational velocities in the HR diagram 

The evolution of v in the HR diagram is shown in Fig. |l3|. 
Starting with models having a velocity of 300 km s _1 on 
the zero-age main sequence, we give some lines of con- 
stant v over the HR diagram. On the MS, we notice in 
particular that the decrease is much faster for the most 
massive stars than for stars with M < 15 M Q . This differ- 
ence remains also present in the domain of B-supergiants. 



During the crossing of the HR diagram, the rotational ve- 
locities decrease fastly, to become very small, i.e. of the 
order of a few km s _1 , in the red supergiant phase. 

It is beyond the scope of this paper to make detailed 
comparisons of the evolution of the distribution of the ve- 
locities over the HR diagram, however, we may notice a 
few points. The fact that the average v is lower for O-type 
stars than for the early B-type stars (Slettebak, 197C ) may 
be the consequence of the higher losses of mass and angu- 
lar momentum in the most massive stars. Also, we remark 
that the increase of v from O-stars to B-stars is larger for 
the stars of luminosity class IV than for class V (Fukuda, 
1982Q. T his is consistent with our models, which show (cf. 
Fig. [13D that the differences of v beween O- and B-type 
stars are much larger at the end of the MS phase. Another 
fact in the observed data is the strong decrease of v for 
the massive supergiants of OB-types. This is predicted by 



all stellar models (cf. also Langer 1998) due to the growth 



of the stellar radii. Further detailed comparisons may per- 
haps provide some new tests and constraints. 
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7. Evolution of the surface abundances 

The chemical abundances offer a very powerful test of in- 
ternal evolution and they give strong evidences in favour 
of some additional mixing processes in massive stars. A 
review of the observations may be found in Maeder and 
Meynet (2000). Here we shall concentrate on the discus- 
sion of the theoretical results and we shall compare them 
with some recent observations. 

The most striking feature appearing in Tables | and | 
as well as on Fig. |l4| is the change of surface abundances 
in rotating stellar models (cf. Langer 1992). The He-, N- 
enrichments and the related C- and O-depletions at the 
surface already occur on the MS. The more massive the 
star is, the more pronounced are the enrichments, sup- 
porting the expectation that mixing becomes more and 



more efficient when the mass increases (Maeder 1998). The 
same is true when the initial rotational velocity increases. 
During the crossing of the Hertzsprung-Russel diagram, 
the evolution is sufficiently short for not allowing any im- 
portant change of the surface abundances (see Fig. |l4| ). 
Further changes occur when the star becomes a red su- 
pergiant and undergoes the first dredge-up. In contrast, 
non-rotating stellar models show no change of the surface 
abundances until the first dredge-up in the red supergiant 
stage (see Table [j]). This implies that these models pre- 
dict some enrichment neither during the MS nor in the 
blue supergiant phase unless a blue loop is formed. More- 
over, as can be seen in Fig. |l4|, the ratios obtained at the 
end of the He-burning phase are significantly lower than 
the ones obtained in rotating models. 



7.1. He-enrichments and v sin i 

Fig. [l5] shows the evolution in the e versus surface veloc- 
ity plane where e = ^h^-h^Hc) > 71 (He) and n(H) being the 
abundances in number of helium and hydrogen at the sur- 
face. The theoretical tracks correspond to the equatorial 
velocities v, while the observed points from Herrero ct al. 
( |1992| , |1999| , |200C| ) are vsini. This implies that the ob- 
served values are smaller on the average by a factor ir/4 
with respect to the theoretical values. 

In Fig. 15, the tracks go generally from the bottom 
on the right to the top on the left. Along a given track, 
the He-enrichment increases when the velocity decreases. 
The shaded zone in Fig. [l5| corresponds to the MS for the 
models with V[ n i — 300 km s _1 . The higher initial mass, the 
greater the He-enrichmcnts which can be reached during 
the MS. We see also that during the MS phase, for initial 
masses inferior to about 60 Mq, all the tracks with v ln i 
= 300 km s _1 follow more or less the same e versus v 



metallicity (Maeder & Meynet 2000b; Meynet 2000a), the 
initial velocity and the age of the star. 

A low surface velocity at a given evolutionary stage 
does not exclude that in the past the star was a fast rota- 
tor. The slow rotation may result from the loss of angular 
momentum by stellar winds and/or from the increase of 
radius of the star. On the other hand, a high velocity in 
the past implies He- and N-enrichments of the surface. 

Let us compare in Fig. [l5| the theoretical tracks in the 
e versus v plane with the observations of OB stars per- 
formed by Herre ro et a l. (|1992|, [1999[ |2000|) . Rece nt works 



(McErlean et al. |l998| ; Smith and Howarth |l998j ) indicate 



that accounting for a microturbulent velocity line broad- 
ening in the model atmosphere reduces the derived He- 
abundances for supergiants later than 09. However, ac- 
cording to Villamariz and Herrero ( 199S ) this effect cannot 
explain all the observed overabundances, especially for the 
earlier types. On Fig. [r] the observed points are plotted in 
the logg c ff versus logT c ff plane where g c ff is the effective 
surface gravity. We estimate g e g = — Q 2 Rsin9 for the 
average orientation angle. Let us note that for the models 
plotted in Fig. [r], there is little difference between the ef- 
fective gravities at the pole and at the equator. Indeed the 
ratio between these two gravities never exceeds 1.3 which 
means a vertical dispersion of about 0.1 dex in Fig. 16. 
Non-rotating and rotating evolutionary tracks are super- 
posed to the observed points in Fig. [r| Since most of the 
enriched stars are in the vicinity of the 120 and 60 M Q 
tracks (see Fig. [l6|), we can wonder wether the changes 
of the surface abundances can be explained as an effect of 
mass loss only. It does not seem to be the case, because the 



relation. However, the relation is changed when the initial 
velocity is different (see the fast rotating 20 and 60 M Q 
models) . The surface He-enrichments on the MS generally 
depend on the following factors: the initial mass, the initial 



part of the track shown on Fig. 16 for the non-rotating 60 
M Q model presents no surface He-enrichment. In the case 
of the non-rotating 120 M Q model only the part of the 
track with logoff inferior to 3.3 has e > 0.14. Thus the 
He-enrichments cannot be accounted for by current evo- 
lutionary models as was already pointed out by Herrero et 
al. ( |1992| , cf. also Maeder 1987). In the following we shall 
suppose that these enhancements are due to rotation. 

Let us consider four groups of stars. In the first group 
we place all the stars presenting no He-enrichment at their 
surface (e < 0.12) and having wsini < 200 km s^ 1 (empty 
squares on Figs, [l^ and 16), in the second one are the 
enriched stars (e > 0.12, solid triangles) having wsini < 
200 km s . The third group (which is empty at present) 
consists of the fast rotators with no He-enrichment. These 
stars would occupy the bottom right corner in Fig. [l5| 
Finally, the fourth group contains the He-enriched stars 
with vs'mi > 200 km s _1 (solid circles). 

Group 1 

The non-enriched stars with low rotation can be in- 
terpreted either as stars with small initial velocities or as 
young fast rotators whose surface has not yet been en- 
ric hed in helium by rotational mixing (see also Herrero et 
al. 1999 ). In this respect let us mention that all the stars 
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Fig. 14. Evolution of the N/C ratios (in mass) at the surface of rotating models for v m i — 300 km s^ 1 . Dotted lines 
joining evolutionary stages with the same value of the N/C ratio are superposed on the MS evolutionary tracks. Values 
of the N/C ratios are indicated at some points during the post-MS phases. The numbers indicated in parentheses 
correspond to the values of the N/C ratio at the surface of the non-rotating models at the end of the He-burning 
phase. The initial ratio (N/C)o is equal to 0.25. For purpose of clarity, only the first part of the 40 M track is shown. 



having log g c ff > 3.7, and which therefore are probably not 
too evolved, present no He-enrichment. 

Group 2 

The most striking feature of the second group of stars 
(e > 0.12 and v s'mi < 200 km s _1 ) is the fact that they are 
distributed in a relatively narrow range of v sin i between 
80 and 160 km s _1 . This may result from the following 
facts: firstly there exists a minimum value of the initial 
velocity for rotational mixing to be able to drive changes 
of the helium surface abundances during the H-burning 
phase. Secondly the observed distribution also reflects the 
way the surface velocity declines during the H-burning 
phase and the narrow range of observed velocities may 
result from some convergence effect as the one mentionned 
in Sect. 6.1 (see also Fig. 12). We see that the high e values 
reached by some of these stars (superior to ~0.16) would 



be compatible with their high initial mass implied from 
their position in Fig. [l6| 

Group 3 

Very interestingly no stars are observed with a v sin i > 
200 km s _1 and no He-enrichment. Keeping in mind that 
the observed stars do not represent a statistical complete 
sample, one can nevertheless wonder why no stars are ob- 
served in this zone. A possibility would be that these very 
fast rotators do not exist, but this is not realistic. Indeed 
stars with velocities as high as 300 km s" 1 and more have 



been observed on the MS (e.g. Penny 1996; Howarth et 
al. 1997; see also Fig. Il5|). Moreover if such stars were not 



formed how to explain the stars in Group 4, having a very 
high v sin i and an important He-enrichment (the solid 
circles in Fig. |l5|) ? These stars are likely the descendants 
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of very fast MS rotating stars (see the discussion in the 
subsection Group 4 below). 

The fact that no stars are observed in the Group 3 may 
indicate that the change of the surface abundances occur 
within a small fraction of the visible MS life. For the fast 
rotating 60 M model plotted on Fig. H], the surface re- 
tains its initial composition (e < 0.11) only during a third 
of its H-burning lifetime. More rapidly rotating models 
would still reduce the fraction of the MS time spent with 
no change of the surface abundances. The lack of fast ro- 
tators with normal surface composition may be due to 
the fact that young massive fast rotators are still embed- 
ded in the cloud from which they formed. The models 
of massive star formation with accretion (Bernasconi and 
Maeder 1996) indicate that when the stars become visi- 



ble, they have already burnt some fraction of their central 
hydrogen, thus some transport of He and N to the stellar 
surface may have occured. 

Group 4 

Could the stars belonging to the fourth group (e > 0.12 
and vsmi > 200 km/ s) be formed by stars previously in 
the low velocity range and which have been accelerated 
for instance by contraction on a blue loop ? The answer is 
likely no. Indeed blue loops cannot accelerate the surface 
beyond the critical velocity which, for a 12 M blue su- 
pergiant at the tip of a blue loop, is of the order of 250 km 
s" 1 . In addition some stars in the group 4 are classified as 
MS stars (see the stars labeled with a V on Figs. |l5| and 
|l6| ). Another possibility would be to consider the stars in 
Group 4 as secondary stars in close binary systems which 
would have been accelerated through the process of mass 
accretion ? But these stars are not observed to belong to 
binary systems. Therefore the most reasonable hypothesis 
is to consider the stars in Group 4 as the natural descen- 
dants of very fast MS rotating stars. 

Their chemical enrichment at the surface is very fast. 
Their chemical structure as a result of the strong rota- 
tional mixing is probably near homogeneity. This view is 
supported by the fact that high values of vsmi are ob- 
served for very high values of e implying that the surface 
velocity does not decrease too much in the course of the 
evolution. This can be accounted for if the star remains 
compact, i.e. in the blue part of the HR diag ram as is the 
case for a strongly mixed star (Maeder 1987 ; see also the 
fast rotating 60 M Q track in Fig. |l5|). Another effect could 
also be important in that respect, i.e. the anisotropy of the 
stellar winds when stars are rotating near break-up. For 
the hot stars, the von Zcipcl ( 1924 ) theorem impli es tha t 
most of the mass is ejected from the pole (Maeder 1999| ). 
This prevents the loss of important angular momentum 
and maintains a high surface velocity. 

If the stars in the fourth group are nearly homoge- 
neous objects, one would expect higher effective gravities 
than observed. Typically, the fast rotating well mixed 60 
Mq model remains at a high value of logoff, at least for 
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Fig. 15. Evolution of the helium surface abundance as 
a function of the equatorial rotational velocity; e = 
ra(H)+rt(He) where n (H c ) and n(H) are the abundances in 
number of helium and hydrogen. All the tracks except two 
were computed for an initial velocity of t>; n i = 300 km s . 
The dotted tracks at the right of the figure correspond to 
fast rotating models, namely a 20 and a 60 M Q model 
with Vi n i equal to 580 and ~400 km s _1 respectively. Only 
the MS of the fast rotating 20 M Q model is shown. The 
observed points (e, usini) are from Herrero et al. ( |1992| , 
1999, 2000| ). The empty squares are for stars with no (or 



very low) He-enrichment at the surface, the solid triangles 
are for objects whose surface present He-enrichments and 
having wsini < 200 km s , the solid circles correspond 
to He-enriched stars with vsmi > 200 km s _1 . A V labels 
the stars of this last group classified as MS stars. 



the portion of the evolution computed here (see Fig. [16]) . 
Moreover, as noted above, the maximum dispersion in 
logffeff due to orientation effects is around 0.1 dex. Could 
the log g c ff be underestimated ? Even if it is difficult on the 
base of the present data to ascertain such a point of view, 
one can note that the sini for these fast rotating stars 
cannot be too far from 1, otherwise one would obtain sur- 
face velocities above the break-up value. This means that 
these stars are seen essentially equator-on. The logg in 
the equator band is inferior to the surface averaged log g 
and thus the observed log g might be underestimated. Of 
course a quantitative analysis requires a detailed study 
of how the variation of log g with the latitude affects the 
spectroscopically determined gravities. 
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Fig. 16. Evolutionary tracks in the logg e g versus logT c ff 
plane where <? e ff is the effective surface gravity. The dashed 
and continuous lines are for the non-rotating and rotating 
(uini = 300 km s _1 ) models respectively. The dotted track 
corresponds to the fast rotating 60 M plotted on Fig. |l5|. 
Only the beginning of the evolution is shown. The posi- 
tions of the stars observed by Herrero et al. (1992, 1999, 



200C ) are indicated with the same symbols as in Fig. |15| 



7.2. Comparisons with the surface abundances of 
supergiants 

Most blue supergi ants present surface enrichments. For 
instance, Walborn ( 1976| , 1988) showed that ordinary OB 
supergiants have He- and N- enrichments as a result of 
CNO processing. Only the small group of peculiar OBC- 
stars has the normal cosmic abundance ratios (cf. also 
Howarth and Prinja 1989; Herrero et al. 1992; Gies and 



Lambert 1992; Lennon 1994). A possibility to explain the 
He- and N-enrichments in supergiants would be that blue 
supergiants are on a blue loop after a first red supergiant 
stage where dredge-up has occurred producing the ob- 
served surface enrichments. However as we shall see below, 
this does not appear as the good explanantion at least for 
some of the observed enrichments. 

Fig. [It] illustrates the changes of the nitrogen to car- 
bon ratios N/C from the ZAMS to the red supergiant 
stage for stars in the mass range from 9 to 20 M©. The 
N/C ratio appears as the most sensitive observable pa- 
rameter. For non-rotating stars, the surface enrichment 
in nitrogen only occurs when the star reaches the red su- 
pergiant phase; there, CNO elements are dredged-up by 
deep convection. The behaviour is the same as for the 
He-enrichmcnts discussed above. For rotating stars, N- 
excesses already occur during the MS phase and they are 



Fig. 17. Evolution as a function of log T e g of Alog^ = 
log(N/C) — log(N/C)i, where N and C are the surface 
abundances (in number) of nitrogen and carbon respec- 
tively, the index i indicates initial values. The continuous 
tracks correspond to rotating stellar models with vi n i = 
300 km s -1 . The dashed and dotted lines are respectively 
for initial Wi n i = 200 and km s . Only parts of the tracks 
with i>ini = 200 km s _1 are shown. The initial masses are 
indicated. The solid squares and triangles correspond to 
observed blue supergiants by Gies and Lambert (1992) 
and Lennon (1994) as reported by Venn (1995b). The solid 
circles show the positions of galactic A-type supergiants 
observed by Venn (1995a, b). 



larger for higher rotation and initial stellar masses. At the 
end of the MS phase of the 12 M model, the N/C ratio 
is enhanced by factors 2.4 and 4.3 for Vi n i = 200 and 300 
km s _1 respectively. These factors are increased up to 4 
and 7.1 for the 20 M models with the same Wi n i- 

On Fig. [I?] we plot also some observations of super- 
giants performed by Gies and Lambert (1992), Lennon 
( |1994| ) and Venn Jl995a| b), The plotted values are 
AlogN/C = [N/C]* - [N/C] B , where the abundance ra- 
tios in number are measured at the surface of the star (*), 
and at the surface of main sequence B-stars (B) supposed 
to have retained their pristine cosmic abundances. Gies 



and Lambert (|99J) provide "LTE" and "NLTE" N/C ra- 
tios. We plot here the smaller "LTE" ratios since, for the 
stars in common with the Lennon sample, they are simi- 



lar to Lennon's results (see the discussion in Venn 1995b) 



From the positions of the observed stars in the log g c s ver- 
sus logT c ff diagram, one obtains that the range of initial 
masses for the A-type supergiants shown on Fig. [l^ are 
between 5 and 20 M Q . 
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From Fig. |l7], one notes first that at the T e g of the 
observed supergiants, non-rotating stellar models predict 
no enrichment at all unless there is a blue loop. However 
it is not likely that all the observed points are at the tip of 
a blue loop. Indeed some of the stars have initial masses 
above 15 M Q and, at solar metallicity, current grids of 



mode ls (Arnctt |1991[ Schall er ct al. |1992| ; Alongi et al. 
1993; Brocato and Castellani 1993| ) only predict blue loops 
for masses equal or lower than 15 M Q . Moreover many of 
the observed stars have N/C ratios too low to result from 
a first dredge-up episode (see for instance the position of 
the blue loop of the non -rotating 9 M Q model in Fig. |i~7|) . 
Therefore Venn ( 1995b| ) suggested that at least those stars 
presenting the lowest N/C ratios are on their way from the 
MS to the red giant branch and have undergone some mix- 
ing in the early stage of their evolution. If such stars are 
not at all accounted for by standard evolutionary tracks, 
rotating models can naturally reproduce their observed 
surface abundances as can be seen on Fig.|l7|. 

Moreover as already noted above, theory predicts 
larger excesses for higher masses, a result in agreement 
with the suggestion of Takeda and Takeda-Hidai (1995) 
recently confirmed by McErlean et al. ( 1999] ). 



8. Conclusion 

Mass loss by stellar winds and rotational mixing in the 
stellar interior are certainly the two hydrodynamical phe- 
nomena which most deeply affect the evolution of massive 
stars. Far from being a small refinement in the physics of 
stellar interior, rotation appears as an essential ingredient 
of future grids of stellar models. In particular among the 
important points wich are not discussed here but which 
will be studied in more details in forthcoming papers are 
the effects of rotation on the population of red and blue 
supergiants at various metallicities, on the evolutionary 
scenarios leading to the formation of Wolf-Rayet stars 
(Maeder 1987; Fliegner and Langer 1995; Meynet 2000b) 
and on the stellar yields (Heger et al. 2000). These ques- 
tions are important for a better understanding of star- 
bursts regions and of the chemical evolution of galaxies. 
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